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ABSTRACT 



Context, It is becoming clear that determination of the abundance of Si using lines of Si n and Si m can lead to quite discordant results 
in mid to late B-type stars. The difference between the Si abundances derived from the two ion states can exceed one dex in some 
cases. 

Aims. We have carried out a study intended to clarify which kinds of B stars exhibit this discrepancy, to try to identify regularities in 
the phenomenon, and to explore possible explanations such as abundance stratification by comparing models to observed spectra. 
Methods. We used spectra from the ESPaDOnS spectropolarimeter and FEROS spectrograph, supplemented with spectra from the 
ESO and ELODIE archives, of magnetic Bp, HgMn, and normal B-type stars ranging in effective temperature from about 10500 to 
15000 K. Using these spectra, we derived abundances using the spectrum synthesis program zeeman which can take into account the 
influence of magnetic fields. For each star, accurate abundances of Si n, Si m, Ti, Cr, and Fe were derived from two separate ~100 A 
windows. Si n abundances were deduced from multiple lines, and Si in abundances were found using AA 4552, 4567, and 4574. 
Results. All magnetic Bp stars in our sample show a discordance between the derived abundances of the first and second ions of 
silicon, with the latter being between 0.6 - 1.7 dex higher. The same behaviour is observed in the non-magnetic stars but to a much 
smaller extent: Si m is enhanced by between 0.3 - 0.8 dex compared to Si n. We do not detect the discrepancy in three stars, HD 22136 
(normal), HD 57608 (HgMn) and HD 27295 (HgMn); these are the only stars in our sample for which the microturbulence parameter 
is significantly different from zero, and which therefore probably have convection occurring in their atmospheres. 
Conclusions. We find that vertical stratification of silicon in the atmospheres of B-type stars may provide an explanation of this 
phenomenon, but our detailed stratification models do not completely explain the discrepancies, which may, in part, be due to non- 
LTE effects. 
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1. Introduction 

The magnetic peculiar A and B-type stars (Ap/Bp) are charac- 
terized by large overabundances of Fe-peak and rare-earth ele- 
ments. They possess strong, non-axisymmetric magnetic fields. 
The observed line-of-sight magnetic field component, (B z ), is 
typically of the order of 1 kG or more in strength, and usually 
varies as the star rotates. The abundance ratios of the chemical 
elements are often variable over the stellar surface ("patchy"), 
leading to line strength and shape variations with stellar rotation. 
This situation is possible because the magnetic field inhibits any 
convective or turbulent mixing. 

In addition to this, [Wade et all d2001bl) argued that most 
magnetic stars exhibit vertical abunda nce variations (chemical 
stratification), an idea first explored by Babel (1992). The phe- 
nomenon was further studied by Rvabchiko va et al.l(l200ll) . who 
compared the chemical abundances derived from the singly- 



* Based in part on observations made with the European Southern 
Observatory (ESO) telescopes under the ESO programmes 067. D- 
0579(A), 072.D-0410(A), 076.B-0055(A), 076.D-0 169(A), 078.D- 
0080(A), and 086.D-0449(A), obtained from the ESO/ST-ECF Science 
Archive Facility. It is also based in part on observations carried out at 
the Canada-France-Hawaii Telescope (CFHT) which is operated by the 
National Research Council of Canada, the Institut National des Science 
de l'Universe of the Centre National de la Recherche Scientifique of 
France and the University of Hawaii. It is based in part on observations 
available on the ELODIE archive. 



and doubly-ionised forms of Pr and Nd in magnetic and non- 
magnetic A-type stars. In pulsating ("roAp") magnetic stars, 
abundances derived from the second ions were found to be 1- 
1.7 dex larger than abundances from the first ions. In contrast, 
they found no significant discrepancies between the abundances 
derived from the first and second ionization states of Pr and Nd 
in the non-magnetic and non-pulsating stars. Lines of Pr m and 
Nd m would normally be formed deeper in the atmospheres of A 
stars than the lines of the singly ionised rare earths, because of 
th e higher local temperature deep in the atmosphere. However, 
as lRvabchikova et all d200ll) pointed out, since electron number 
density decreases outward, the second ionization state can also 
dominate in the very uppermost atmospheric layers. Because of 
this effect, they were able to understand the discrepancy between 
ion states as being due to large rare earth overabundances con- 
centrated in a thin layer high in the atmosphere. 

A related discrepancy has be en observed for abund ances of 
Si derived from two ion states. ISemenko et al.l (120081) studied 
the magnetic Bp star HD 45583, which has effective tempera- 
ture T e g = 13000 K, \ogg = 4.0 and a strong magnetic field, 
with (B z ) varying from about -3500 to +4000 G. Different abun- 
dances must be assumed for Si n and Si m to adequately fit the 
theoretical to the observed line profiles. They interpret this dis- 
crepancy as being due to vertical stratification of silicon in the 
atmosphere, but point out that increasing T e g from the assumed 
value also brings the abundances derived from lines of Si n and 
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Si m into better agreement. Bail ey et all (1201 2b found a simi- 
lar discordance in the derived abundances of the first and sec- 
ond ionization states of Si in HD 133880, another strongly mag- 
netic Bp star «fl z ) varies from about -4000 to +2000 G) having 
Teff = 13000 K and logg = 4.34. They find that an enhancement 
of about 1 dex above the Si n abundance is necessary to ade- 
quately fit the Si m lines. Since non-LTE effects in the silicon 
lines are expected to be small below a T e g of about 15000 K, 
they suggest that the most probable cause is strong stratification 
of Si in the atmosphere. 

The study of vertical abundance stratification in magnetic Ap 
stars has considerably expanded in the past fe w years with work 
led b y Koc h ukhov and Ryabchikova (e.g. Kpchukhov et all 



Table 1. Stars analyzed in this study. Listed are the star desig- 
nations, instrument used, spectral resolution, spectral range and 
the S/N. 



2011). This 



[20061 120091 IShulvaketall 120091 iPandev et all 
group has developed refined methods of recovering the vertical 
distribution of several elements in a stellar atmosphere by simul- 
taneous fits to a number of spectral lines of different strengths, 
excitation potentials, and ionisation states. They have studied 
several cool magnetic Ap stars in detail, and mapped the verti- 
cal distribution of a number of elements in each, in several cases 
including Si. Typically, it is found that Si has high abundance 
deep in the atmosphere, and an abundance that declines strongly 
towards low optical depths, sometimes by as much as 3 dex over- 
all. 

Up to now, the study of discrepancies among different ion 
states, and of stratification as a probable explanation, has mainly 
been confined to cool magnetic Ap stars. It is time for a first 
survey of higher effective temperatures. In this paper, we report 
an empirical study of abundances derived from lines of Si 11 and 
Si in in various classes of B-type stars including magnetic Bp, 
HgMn, and normal B stars. Stars were selected from the ESO, 
CFHT, and ELODIE archives (as well as from our own data from 
previous CFHT and ESO programmes). The stars studied have 
r e ff values that range from about 10500 to 15000 K, and have 
low enough projected rotational velocity v sin i (less than about 
100 km s _1 ) that lines of Si in can be unambiguously detected. 
The following section discusses the spectroscopic observations 
used. Sect. 3 outlines the modelling technique and the spectrum 
synthesis program used. Sect. 4 presents the abundances we de- 
rive from lines of Si 11 and Si m. Sect. 5 investigates the possible 
effects of non-LTE in deriving abundances of Si 11 and Si m with 
increasing T e ff and we test whether the observed discrepancies 
can be accounted for using stratified abundance models. Sect. 6 
summarises and discusses the work we present. 

2. Stellar sample 

We wish to study the extent to which the abundance of Si as de- 
rived from Si 11 and Si m is concordant or discordant. With the 
tools readily available to us, this can only be done in a limited 
range of r e ff values or spectral range. Our low-temperature limit 
is set by the fact that no Si in lines are visible in optical spectra 
for r e ff less than about 1 1 000 K. The high-temperature limit is 
set by the gradual ons et of non-LTE effects at about 15 000 K 
dPrzvbillaet all 1201 lh . This upper limit to our study, which is 
not well-defined, is imposed by two main factors. First, the non- 
LTE codes developed up to now do not incorporate radiative 
transfer in the presence of a magnetic field. Thus, they cannot 
be used directly to study magnetic stars which are an important 
class of objects for this study. Furthermore, as we will show be- 
low, the publis hed grid of non-LTE eq uivalent widths for Si 11 
and in lines of iBecker & Butlerl (Il990al) do not even approxi- 
mately reproduce the line ratios in the normal B3 IV star 1 Her 
= HD 160762, and recent non-LTE abundance studies of early 



Star 


Instrument 


R 


A ( A) 


S/N 




normal stars 






HD 22136 


ESPaDOnS 


65000 


3690-10481 


504 


HD 160762 


ESPaDOnS 


65000 


3690-10481 


600 


HD 162586 


UVES 


110000 


3070-10398 


655 


HD 170054 


ESPaDOnS 


65000 


3690-10481 


283 


HD 179761 


ELODIE 


42000 


4000-6800 


125 


HD 195810 


ESPaDOnS 


65000 


3690-10481 


774 


HD 222173 


ELODIE 


42000 


4000-6800 


159 




HgMn stars 






HD 27295 


UVES 


110000 


3750-6800 


218 


HD 57608 


FEROS 


48000 


3528-9217 


142 


HD 78316 


UVES 


110000 


3750-6800 


265 


HD 175640 


UVES 


110000 


3750-5800 


439 


HD 178065 


UVES 


110000 


3750-5800 


461 


HD 186122 


ESPaDOnS 


65000 


3690-10481 


662 


HD 193452 


CF4 


120000 


4630-5478 


556 




magnetic Bp stars 






HD 10840 


UVES 


110000 


3070-10398 


824 


HD 45583 


FEROS 


48000 


3528-9217 


231 


HD 47116 


UVES 


110000 


3070-10398 


527 


HD 49333 


UVES 


110000 


3070-10398 


791 


HD 61045 


ESPaDOnS 


65000 


3690-10481 


298 


HD 74168 


UVES 


110000 


3070-10398 


334 


HD 74535 


FEROS 


48000 


3528-9217 


310 


HD 133880 


ESPaDOnS 


65000 


3690-10481 


365 


HD 137509 


UVES 


110000 


3690-10481 


430 


HD 147010 


ESPaDOnS 


65000 


3690-10481 


587 


HD 199728 


UVES 


110000 


3070-10398 


593 


HD 223640 


UVES 


110000 


3070-10398 


560 


HD 304842 


FEROS 


48000 


3528-9217 


262 


HD 318107 


ESPaDOnS 


65000 


3690-10481 


320 


BD+00 1659 


ESPaDOnS 


65000 


3690-10481 


475 


BD-19 5044L 


ESPaDOnS 


65000 


3690-10481 


358 


BD+49 3789 


ESPaDOnS 


65000 


3690-10481 


191 



B stars appear to systematically reject use of the Si 11 lines. It is 
not at all clear that current non-LTE codes are able to correctly 
compute line ratios for Si 11 lines. Therefore we try to limit our 
study to effective temperatures low enough that strong non-LTE 
effects are not expected; this limiti ng temperature is conven tion- 
ally set at r eff = 15000 K (e.g. iNieva & Przvbillal [20F2T) . We 
have included one hotter star, HD 160762, in order to address 
non-LTE questions later in the paper. 

The spectra for this study were obtained from a variety of 
sources. A large portion of the stars studied were acquired from 
previous observing runs using the ESPaDOnS spectropolarime- 
ter (with resolving power R = 65000) at the Canada-France- 
Hawaii Telescope (CFHT), and the FEROS spectrograph (R = 
48000) at the European Southern Observatory's (ESO) La Silla 
Observatory. The remaining spectra used were acquired from the 
ESO, CFHT and ELODIE (Observatoire de Haute Provence) 
archives. We required that each individual spectrum have a 
signal-to-noise ratio (S/N) greater than about 100 and spectral 
resolution no less than about 40 000. As discussed above, we 
selected stars of between about 10500 and 15000 K, requir- 
ing also that they have sufficiently slow rotation (v sin i less than 
about 100 km s _1 ) that at least some individual lines are not 
strongly blended. 

In the effective temperature range of interest for this study, 
there are (at least) three major classes of main sequence stars. 
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Fig. 1. Blackwell diagrams (from left to right) for HD 57608, 
HD 186122, HD 179761 and HD 222173. Shown are lines for 
Fe ii AA 4541, 4555, 4576, 4583, 5018 and sometimes 6120 
and/or 6129. Each panel spans 1.0 dex. 



In addition to normal late B stars (which generally do not have 
detectable magnetic fields, and have approximately solar sur- 
face abundances), we also find magnetic Bp stars, which usually 
show quite marked abundance differences compared to the Sun, 
and the non-magnetic HgMn stars, which are typically found 
in fairly close binary systems and show quite non-solar over- 
abundances of several low-abundance elements such as P, Mn, 
Ga, and Hg. We have included several examples of each of these 
classes in our sample. The sample of stars utilized in this study 
is summarised in Table Q] 

The fundamental parameters T e ff and logg were derived 
for each star using Geneva and uvby/3 photometry. For the 
Geneva photometry , we use the Fortran program described by 
iKunzli et alJ dl997h . For stars with Stromgren uvby/3 photometry, 
the Fo rtran program "UVBYBETANEW" of iNapiwotzki et aH 
d!993l) . Our version of this program corrects the T e g of the 
magnetic Ap/Bp stars to a suitable Ap temperature scale (see 
lLandstreet et all 120071) . For stars for which both sets of pho- 
tometry were available, the average value was ta ken. We were 
able t o compare our T e fi values with those of Netopil et al. 
(2008) in a few cases; our values agree with theirs except for 
HD 13388 0, for which T e s is somewhat more uncertain, as dis- 
cussed by iBailev et~aTI (120121) . Finally, we have estimated the 
u ncertainty in T e ff t o be about 500 K, following the discussion 
of lLandstreet et all d2007l) . The uncertainty in log g is estimated 
from the level of agreement between values derived from Geneva 
and Stromgren photometry to be about 0.2 dex. The physical 
param eters for HD 160762 were taken from iNieva & Przvbillal 
(2012), and have uncertainties that are about half as large as 
those of the other normal stars. Table ^summarises all the phys- 
ical parameters of the stars modelled in this study, including pa- 
rameters discussed later in this paper. 



3. Modelling spectra 

3.1. Modelling with the spectrum synthesis code zeeman 

Modelling and abundance determination was carried out using 
the Fortran spectrum synthesis code zeeman (se e lLandstreet, 
1988; lLandstreet etail Il989t IWadeetall 12001 at IBailev et all 
2012). zeeman accepts as input an assumed magnetic field geom- 
etry for stars with magnetic fields (such as the magnetic Ap/Bp 
stars), but functions as a conventional li ne synthesis program fo r 
modelling of non-magnetic stars (e.g. lLandstreet et al.L [2009). 
The program assumes LTE atmospheric structure and line for- 
mation. An appropriate stellar atmosphere model is interpolated 
in a grid of ATLAS 9 solar abundance models. For this study, a 
uniform abundance is assumed over the stellar atmosphere even 
for stars such as magnetic Bp stars for which this may be only a 
rough first approximation. 

zeeman functions by comparing a synthesized spectrum to an 
observed one in one or more selected wavelength windows. By 
iteratively minimizing the mean square difference between the 
computed and observed spectrum, zeeman automatically deter- 
mines the best fit radial velocity vr and the best value of v sin i 
for the spectral windows modelled. It optimises the fit to all the 
lines of a single chemical element within the selected window(s), 
determining the value of the abundance of that element relative 
the H. This process is repeated for other elements until most of 
the lines in each spectral window studied are modelled as well 
as possible with the limitations imposed by the assumptions built 
into the model. The adopted abundances are obtained from the 
average of two ~100 A windows, with the uncertainty estimated 
from the observed scatter of the computed values in the modelled 
spectral windows. 

For non-magnetic stars, the microturbule nce paramet e r £, wa s 
determined using the method described by lLandstreetl (Q998 ). 
The abundance is iteratively deduced for a number of individual 
lines of Fe for a small grid of assumed £ values. The resulting 
abundance values are plotted in a Blackwell diagram for each 
line as a function of The value of £ for which the abundance 
values from the lines studied agree best is accepted as the op- 
timal value, and the corresponding value of Fe/H abundance is 
adopted. Examination of the dispersion of the abundance values 
as a function of £ yields estimates of the uncertainties of both Fe 
abundance and Several examples of this procedure are shown 
in Figure [1] 

In most cases, the microturbulence parameter does not seem 
to be significantly different from zero, but we think that we de- 
tect a non-zero value in three stars: the normal star HD 22136, 
and the two HgMn stars HD 57608 and HD 27295. All three 
of these stars are near the low temperature limit of our sample, 
and all three have rather sharp lines, making abundance analysis 
more precise and perhaps making detection of non-zero £ easier. 
Detection of non-zero microturbulence is quite interesting, as it 
suggests the presenc e of convective mixing in at least part o f 
the atmosphere (e.g. lLandstreetl 1 19981; lLandstreet et al.L |2009). 
This would in turn suggest that the atmosphere is at least partly 
mixed, so that a stratified abundance distribution could be de- 
stroyed by rapid mixing, if the mixing extends through enough 
of the atmosphere. 

For magnetic stars, the microturbulence parameter £ was set 
to 0, as the magnetic field is thought to suppress convective mo- 
tions. 

For the magnetic stars HD 137509, HD 318107, and 
HD 133880, detailed colinear multipole magnetic field ge- 
ometries have been reported in the literature that could 
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Table 2. Physical parameters for the stars in this study. We adopt a uniform uncertainty of ±500 K in T e s and +0.2 in log g. 



Star 


r cff (K) 


logg 


v sin i (km s 1 ) 


* O p O 


B rf (G) 


B„ (G) 


B D « (G) 


£(kms-') 








normal stars 












HD 222173 


11800 


3.4 


67 + 5 






- 


- 


- 


< 1.4 


HD 22136 


12700 


4.2 


15 ±2 






- 


- 


- 


1.1 ±0.3 


HD 162586 


12700 


4.0 


28 ±2 






- 


- 


- 


< 1.6 


HD 179761 


12900 


3.4 


14 ±2 






- 


- 


- 


< 1.5 


HD 195810 


13700 


3.7 


52 ±4 






- 


- 


- 


< 1.1 


HD 170054 


14500 


4.3 


25 ±2 




_ 


- 


- 


- 


< 1.2 


HD 160762 


17500 


3.8 


6± 1 






- 


- 


- 


1 ± 1 








HgMn stars 












HD 193452 


10600 


4.1 


1.4 ± 1 






- 


- 


- 


< 1.0 


HD 57608 


10900 


3.1 


5.0 ± 1 






- 


- 


- 


1.3 ± 0.4 


HD 27295 


11800 


4.2 


4.9 ± 0.3 






- 


- 


- 


1.8 ±0.5 


HD 175640 


12000 


4.0 


1.5 ± 1 






- 


- 


- 


< 1.7 


HD 178065 


12300 


3.6 


2.0 ± 1 






- 


- 


- 


< 1.5 


HD 186122 


12900 


3.7 


1.0 + 0.5 






- 


- 


- 


< 1.4 


HD 78316 


13400 


3.9 


6.8 + 0.5 






- 


- 


- 


< 1.6 








mag 


netic Bp stars 










HD 47116 


11000 


4.1 


30 + 2 






500 


- 


- 





HD 10840 


11600 


3.6 


35 + 5 






500 


- 


- 





HD 318107 


11800 


4.2 


7 ± 1 


22 


65 


25600 


-12800 


900 





HD 199728 


12200 


3.7 


62 + 6 






800 


- 


- 





HD 223640 


12300 


4.4 


31 + 3 






500 


- 


- 





BD+00 1659 


12500 


4.0 


7.0 ± 1 






1200 


- 


- 





HD 304842 


12500 


3.9 


65 + 5 






100 


- 


- 





HD 45583 


12700 


4.2 


70 + 6 






8000 


- 


- 





BD-19 5044L 


12800 


4.5 


15 + 3 






800 


- 


- 





BD+49 3789 


12900 


4.2 


85 + 5 






1700 









HD 74168 


12900 


4.5 


63 + 4 






200 









HD 61045 


13000 


4.1 


64 + 3 






1300 









HD 147010 


13000 


4.4 


15 + 2 






15000 









HD 133880 


13000 


4.3 


103 ± 10 


55 


78 


-9600 


-23000 


1900 





HD 137509 


13100 


4.3 


20 + 2 


81 


64 


3100 


41900 


-500 





HD 74535 


13600 


4.3 


45 + 4 






300 









HD 49333 


15100 


3.9 


68 + 3 






800 










be used for modelling: (see iKochukhovl l2006t iBailev et all 
1201 ll 1201 2l respectively). For other stars, the line-of-sigh t 
magnetic field measurements, ( B 7 ), of |Bagnulo et alj (12 006): 
iKochukhov & Bagnuld d2006h : lLandstreet et alJ d2008t) and 
available unpublished data were used to estimate the polar field 
strength of a simple dipole field model. In these cases, a dipolar 
field of about three times the root-mean square of the available 
(B z ) measurements was used in computing the synthetic spectra. 
For stars for which abundances of Fe-peak and rare earth ele- 
ments clearly indicate a Bp nature, but for which no (B z ) mea- 
surements are available, a generic polar field of ~500 G was 
used. For all stars for which no previous magnetic modelling 
was available, we set the inclination of the rotation axis to the 
line of sight, i, and the obliquity of the magnetic field axis to the 
rotation axis, ft, to zero (z'=/J=0). 

3.2. Atomic data 

The atomic data (energy levels, oscillator strengths, damp- 
ing constants, and Lande factors) used for this analysis 
were taken from the Vienna Atomic Lin e Database (VALD) 
dPiskunov et all 11995b iRvabchikova et all 1 19971 iKupkaet al. . 
[19990 2000). When Lande factors were not available from this 
source, they were computed assuming LS coupling. 

Since it is clear that this study depends sensitively on the 
accuracy of the atomic line data, particularly oscillator strengths, 
we have endeavoured to test the VALD gf values extensively. 



One test of the Si n gf values is to use them for determina- 
tion of the abundance of silicon in the star Sirius = HD 48915, 
for which very accurate fundamental and other parameters (e.g. 
the value of ^) are available, and for which a number of excel- 
lent UVE S high-resolutio n spectra are available from the ESO 
archive (Lands treetll201 ll) . Abundances of Si n were determined 
from the lines AA 4621, 5041 & 5056, 5957 & 5978, and 6347 & 
6371. The derived abundances have a mean of log A^si / log A^h = 
-4.25 ± 0.07 dex, where the uncertainty is the dispersion of the 
four values. The value, about 0.2 dex larger than the solar abun- 
dance, is not unexpected in this hot Am star, and the dispersion 
indicates that the gf values used are closely compatible. 

This test has also been carried out on several of the stars in 
our sample, to see how consistently the VALD data give the same 
abundances with different spectral lines and line groups. To ob- 
tain a realistic idea of the uncertainties in deriving abundances 
for Si ii, we select a sample of magnetic, normal, and HgMn stars 
that are well spread in T e ff (within the limitations of the data at 
hand). For these stars we derive abundances from /14621 individ- 
ually, /L15041, 5055 and 5056 simultaneously, and AA5957 and 
5978 simultaneously. Table [3] presents our results, line by line, 
and with the computed dispersion of the three abundance values 
for each star. In general, the abundances we derive from the var- 
ious lines of Si n show dispersion which is only slightly greater 
than that found for Sirius. The magnetic stars show the greatest 
tendency for the dispersion of abundances derived from different 
lines of Si n to be larger than we found for Sirius; however, this 
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Table 3. Derived abundances from individual lines of Si n for a selection of stars in our sample. 



HgMn stars Normal stars Magnetic stars 

HD HD HD 



Line (A) 


178065 


27295 


57608 


78316 


222173 


22136 


170054 


162586 


74168 


49333 


223640 


47116 


4621 

5041-55-56 
5957-78 


-4.67 
-4.56 


-4.52 
-4.58 
-4.54 


-4.89 
-4.78 
-4.66 


-4.34 
-4.44 
-4.26 


-4.66 
-4.63 
-4.48 


-4.67 
-4.64 
-4.59 


-4.64 
-4.81 
-4.70 


-4.67 
-4.82 
-4.72 


-4.00 
-4.11 

-3.93 


-3.90 
-4.15 
-3.91 


-3.09 
-3.20 
-2.98 


-3.21 
-3.48 
-3.27 


Dispersion 




0.03 


0.12 


0.09 


0.10 


0.04 


0.09 


0.08 


0.09 


0.14 


0.11 


0.14 



result is expected, as the magnetic stars typically have horizon- 
tally non-uniform abundances and may also have vertically strat- 
ified abundance distributions (Wade et al., 2001b). Nevertheless, 
we note that in general assuming a uniform abundance of Si fits 
all lines of Si n in all stars analyzed within a fitting uncertainty 
of about +0. 15 dex or better. 

A further test of the atomic data is to compare the VALD 
gf values with the data i n the comprehensive NIST co mpilation 
of atomic data for Si bv lKelleher & Podobedoval d2008h . In this 
data set, the recommended atomic data are selected on a line-by- 
line basis from available experiments and computations, without 
consideration of internal homogeneity. Most of the log gf val- 
ues of Si n needed for our work are estimated by those authors to 
have uncertainties of around +0. 1 dex. Overall, for the lines stud- 
ied by us, there is only a very small difference (NIST - VALD) 
in the mean of all required log gf values of about -0.04 dex 
for Si n, and +0.1 dex for Si m. As we shall see below, this is 
not a large enough difference to play an important role in mea- 
surements of the difference in abundances found in some stars 
between Si n and Si m, although if we used the NIST gf values 
instead of VALD data, this would tend to decrease the abundance 
differences found by ~ 0.15 dex. 

We have repeated the Si abundance determinations for Sirius 
with the NIST Si n data. The average abundance found is in this 
way is -4.20 + 0.14. The change in the mean is thus very small, 
but the dispersion from line group to line group is doubled with 
respect to the results found using the VALD data. This kind of 
difference has been found in earlier comparisons of NIST and 
other, more homogenized, data (e.g. lSigut & Landstreeti ll990). 
and appears to arise from the NIST preference for data selection 
for individual lines from single data sources. In particular, the 
NIST gf value for the Si n line at 5041 A is quite discrepant 
with the VALD value, and when the NIST values are used for 
synthesis of the spectral region containing the AA 5041 and 5056 
lines, the best fits are clearly not concordant between the two 
lines. It is because of this larger dispersion that we prefer to use 
the VALD atomic data. 

For Si in we are able to determine the Si abundance using 
only the lines of multiplet (2) at 4552, 4567, and 4574 A. The 
log gf values adopted by VALD, and by us, for these three lines 
are all about 0. 1 dex smaller than those in the NIST compila- 
tion. The relative strengths of the lines are closely identical in 
both datasets. For consistency with our Si n analyses, we use the 
VALD values. 

We have also double-checked that our assumed ioni- 
zation potentials are very c losely the same as those of 
iKelleher & Podobedoval (120081) . In summary, there is no obvi- 
ous reason to think that the results we find below are due to some 
major defect in the available spectral line data. 



3.3. Si abundance uncertainties resulting from errors in 
stellar parameters 

The uncertainty in log g and especially in T e ff introduces signif- 
icant uncertainty in the derived abundances, especially of Si m. 
This issue has been explored for the different classes of stars in 
this study by repeating the determination of the best fitting abun- 
dance of several stars with values of T e ff or log g altered by the 
adopted uncertainties, and comparing the resulting abundance 
values with those derived with the adopted stellar parameters. 

Over most of the temperature range studied, Si n is the dom- 
inant form of Si throughout much of the atmosphere. Therefore, 
changing the value of T e ff by 500 K typically results in a fairly 
small change in the abundance of Si n which best fits the ob- 
served lines. Generally, the change in derived abundance is of 
the order of 0.05 dex, with the derived abundance decreasing 
slightly if the adopted T e g value is increased. Apparently the in- 
creased population of the high excitation lower levels of the Si n 
lines (at around 10 eV) more than compensates for the effect of 
the Saha equation in decreasing the fraction of Si in the form of 
Si n. The changes to the abundance derived from lines of Si n for 
a change in log g of 0.2 dex is also of order 0.05 dex, with an in- 
crease in log g resulting in a decrease in the abundance deduced 
from Si n, as the Saha equation forces more Si from the doubly 
into the singly ionised state. 

In the temperature range studied, the fraction of Si ra in- 
creases strongly from a very small fraction at the low tempera- 
ture end to about half of the total Si near 15000 K. Furthermore, 
the lower levels of the lines of multiplet (2) are at almost 19 eV. 
As a result, the abundance of Si derived from lines of Si in is a 
much stronger function of T e g . With a change of T e g of +500 K, 
the deduced Si abundance can change by as much as roughly 
-0.3 dex, decreasing somewhat towards the high end of the tem- 
perature range studied. The deduced abundance also changes 
more with a change in log g of 0.2 dex than is the case for Si n; 
the derived abundance can change by an amount of order 0. 1 dex . 



4. Measured Si n and m abundances in the B star 
sample 

In this Section we tabulate and discuss the results of our abun- 
dance analysis of the stars of the sample discussed above. 

Table |4]presents the final mean abundances we find for each 
star. The first four columns list the star name (generally an HD 
Number), and recall the T e g, v sin i, and £ values from Table [2] 
Next we tabulate abundance values found for Ti, Cr, and Fe. The 
abundances of these three elements had to be determined in or- 
der to evaluate and test the microturbulence parameter indepen- 
dently of the Si lines (for the non-magnetic stars only), to cor- 
rectly include the effects of weak blends with the Si lines mod- 
elled, and to allow a precise determination of vr and v sin i. For 
all stars, the abundance of iron was found from simultaneously 
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Table 4. Abundances derived for Ti, Cr, Fe, and the first and second ionization states of Si. Abundance uncertainties are estimated 
to be about +0.2 dex for the magnetic Bp stars and +0.1 dex for the non-magnetic normal and HgMn stars. Exceptional values 
are indicated in parentheses. The second to last column indicates the difference between the Si n and Si in abundances and the last 
column denotes the significance ASi/<x D of the difference. 



Star 


T sB (K) 


v sin i (km s ') 


f (kms-'j 


log(Ti/H) 


log(Cr/H) 


log(Fe/H) 


log(Si/H) 


ASi 


0~D 


















Si ii 


















normal stars 












HD 222173 


11800 


67 ±5 




< 1.4 


-7.80 


-6.64 


-4.96 


—A SQlTI 1 


—A fnCfl 7 s ! 


0.56 


2.2 


HD 22136 


12700 


15 + 2 


1.1 


±0.3 


-7.32 


-6.34 


-4.82 


—4 f.1 


-4.41 


0.22 


1.4 


HD 162586 


12700 


28 ±2 




< 1.6 


-6.95 


-6.26 


-4.84 


-4.74 


-3.90 


0.84 


5.9 


HD 179761 


12900 


14 ±2 




< 1.5 


-7.43 


-6.38 


-4.75 


-4.23 


-4.72 


0.49 


3.5 


HD 195810 


13700 


52 ±4 




< 1.1 


-7.41 


-6.35 


-4.58 


-4.24 


-4.65 


0.41 


2.9 


HD 170054 


14500 


25 ± 2 




< 1.2 


-7.46 


-6.67 


-4.95 


-4.72 


-4.21 


0.51 


3.6 


HD 160762 


17500 


6± 1 




1 ± 1 






-4.49 


-5.34 


-4.10 


1.24 


8.9 












Hg 


Mn stars 












HD 193452 


10600 


1.4 ±1 




< 1.0 


-6.50 


-5.73 


-4.22 


—A 70 








HD 57608 


10900 


5.0 ± 1 


1.3+0.4 


-8.08 


-6.81 


-5.31 


—A 78 


—A 87 


-0.04 


0.28 


HD 27295 


11800 


4.9 ± 0.3 


1.8 


±0.5 


-5.83 


-6.03 


-5.31 


—A S8 


—A 4.1 


0.17 


1.2 


HD 175640 


12000 


1.5 ± 1 




< 1.7 


-6.54 


-5.79 


-4.83 


—A Af, 


— 1 If, 

— J. / U 


0.70 


5.0 


HD 178065 


12300 


2.0 ± 1 




< 1.5 


-6.49 


-6.04 


-4.65 


—A fi1 


—A 7Q 


0.32 


2.3 


HD 186122 


12900 


1.0 ±0.5 




< 1.4 


-6.68 


-7.38 


-4.09 


-"i SA 


—A QO 


0.64 


4.5 


HD 78316 


13400 


6.8 ± 0.5 




< 1.6 


-6.79 


-6.06 


-4.45 


—A 


— ^ 01 

J.71 


0.44 


3.2 












magnetic Bp stars 












HD 47116 


11000 


30 ± 2 







-4.65 


-4.93 


-3.65 


1 ^7 

J.jZ 


7 

Z. jy 


0.93 


3.3 


HD 10840 


11600 


35 ± 5 







-6.68 


-5.56 


-3.49 




7 71 
— Z. / 1 


0.77 


2.7 


HD 318107 


11800 


7 ± 1 







-5.05(0.1) 


-4.50(0.15) 


-3.00(0.2) 




7 I'X 
— Z. / J 


0.92 


3.7 


HD 199728 


12200 


62 ± 6 







-7.38 


-5.33 


-4.07 


J.U7 


7 Af\(C\ 
Z.^fu^U.Zj ) 


0.63 


2.0 


HD 223640 


12300 


31 ±3 







-6.89 


-4.65 


-3.42 


-3.42 


-2.17 


1.25 


4.4 


BD+00 1659 


12500 


7.0 ± 1 







-6.69 


-5.02 


-3.79 


-3.53 


-2.51 


1.02 


3.6 


HD 304842 


12500 


65 ±5 







-6.82 


-6.16 


-4.34 


-3.65 


-2.91(0.25) 


0.74 


2.3 


HD 45583 


12700 


70 ±6 







-5.67 


-4.68 


-3.40 


-3.33 


-2.10(0.25) 


1.23 


3.8 


BD-19 5044L 


12800 


15 + 3 







-6.70 


-5.80 


-4.15 


-3.99 


-2.76 


1.23 


4.3 


BD+49 3789 


12900 


85 ±5 







-6.23 


-5.48 


-4.09 


-3.44 


-2.59(0.3) 


0.85 


2.4 


HD 74168 


12900 


63 ±4 







-6.77 


-6.30 


-4.81 


-4.01 


-2.95(0.25) 


1.06 


3.3 


HD 61045 


13000 


64 ±3 







-5.96 


-5.24 


-3.90 


-3.79 


-2.56(0.25) 


1.23 


3.8 


HD 147010 


13000 


15 + 2 







-5.81 


-4.18 


-3.46 


-4.03 


-3.46 


0.57 


2.0 


HD 133880 


13000 


103 ± 10 







-4.86 


-4.49 


-3.30 


-2.68 


-1.23(0.3) 


1.45 


4.0 


HD 137509 


13100 


20 ± 2 







-4.62 


-4.20 


-3.19 


-3.52 


-2.02 


1.50 


5.3 


HD 74535 


13600 


45 ±4 







-6.25 


-5.60 


-3.98 


-4.25 


-2.53 


1.72 


6.1 


HD 49333 


15100 


68 ±3 







-5.97 


-5.66 


-4.18 


-3.99 


-2.50(0.3) 


1.49 


4.1 



fitting three lines of Fe n: 4541, 4555, and 4583 A. Similarly, 
where possible we derived the abundance of chromium from 
lines of Cr n at 4558, 4588, and 4592 A and titanium from lines 
of Ti ii at 4565 and 4572 A. The abun dance of Fe and the valu e 
of £ for HD 160762 were taken from lNieva & PrzvbiHal(l2012l) . 
but the Si abundances were determined by us in the usual way. 

For the magnetic Bp stars, we estimate the uncertainties of 
the abundances of the Fe-peak elements to be ±0.2 dex, a typical 
deviation from the derived abundance necessary to produce an 
unsatisfactory fit of the model to the observed spectra. The un- 
certainty derived in the same manner for the non-magnetic stars 
suggests a smaller uncertainty of ±0.1 dex. Uncertainties that 
differ from the aforementioned values are noted in parentheses. 
This is the case for stars with v sin i greater than about 60 km s _I 
where line blending makes the derived abundances somewhat 
less certain. 

Following the Fe-peak elements, two columns list abun- 
dances derived from lines of Si n and Si m separately, for each 
star in our sample. The next column gives the difference of the 
two Si abundance determinations, in the sense Si m - Si n, and 
the final column gives the difference divided by the estimated 
uncertainty of the difference, including only the fitting uncer- 
tainties discussed in this section. (Recall that the uncertainty in 



T e ff substantially increases the uncertainty in the abundance de- 
rived from Si m, to at least ±0.3 dex.) 

It will immediately be seen that the values of \og{Ns, JNu) 
are quite discordant as derived from the two different ionisation 
states of Si. The reader may wonder how robust these large dif- 
ferences are. Figure [2] illustrates the type of fits we obtain with 
magnetic and non-magnetic stars using abundances derived from 
lines of Si n and Si m. The situation for the lines of these two 
ions is mostly quite different than the fits found for the Fe-peak 
elements Ti, Cr and Fe, for which all lines of the elements are 
generally fit well with a single value of the abundance. In both 
the HgMn star (upper spectrum) and the magnetic star (lower 
spectrum), the abundance that best fits the lines of multiplet (2) 
of Si in at 4552-67-74 A (model fit in blue) predicts a Si n line 
at 4621 A that is far too strong. Conversely, the abundance that 
fits the Si ii line at 4621 A, and the other Si n lines modelled, 
leads to model lines for the Si m lines that are much too weak. 
This discrepancy is really striking for the magnetic Bp star. It is 
present but at a less spectacular level for most of the HgMn stars, 
and even seems to be present for some of the normal stars. 

The sharp-lined HgMn stars, which have a roughly solar 
abundance of Si, allow us to establish the lower limit of 
for which Si m is observable. This was done by searching for 
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4552 4553 4567 4568 4569 4574 4575 4620 4622 4624 
Wavelength (A) 



Fig. 2. Model fits to lines of Si n (right) and Si in (left three 
panels) to the non-magnetic HgMn Star HD 78316 (upper spec- 
trum) and the magnetic Bp star BD-19 5044L (lower spectrum). 
The black lines are the observed spectra, red lines are the best fit 
model using the abundance derived from Si n and the blue lines 
are the best fit model using the abundance derived from Si m. 



the presence of /14552, the strongest Si m line in multiplet 2, 
which also includes AA4561, and 4574. We are able to detect a 
very weak line of Si in in HD 57608 (T ef{ = 10900 K), however, 
in a slightly cooler star, HD 192452 (T eS = 10600 K) we are 
unable to unambiguously detect any lines of Si m. This thresh- 
old also depends somewhat on log g, but for a solar abundance 
of Si, we can begin to detect Si in in main sequence stars at 
r eff * 11000K. 

5. Possible explanations of the Si n-ra anomaly 

The results presented in the previous section show that the abun- 
dances derived from the singly and doubly ionised states of Si 
are in general different. This difference is always of the order 
of one dex for the magnetic stars in our sample, and appar- 
ently varies from near zero to almost one dex for the normal 
and HgMn stars. The difference is always in the same sense: if 
the abundance is chosen to have a value consistent with the Si n 
lines, the predicted strengths of the Si m lines are always too 
small. 

It is first worth noticing that there are three stars for which no 
hint of discrepancy between the abundance of Si derived from 
the two ionization states. These are exactly the three stars for 
which we believe that we detect a non-zero microturbulence pa- 
rameter and thus presumably the presence of atmospheric con- 
vection. It is not obvious what makes these stars different from 
other normal or HgMn stars of similar fundamental parameters, 
since we expect the presence or absence of atmospheric convec- 
tive zones to vary smoothly with T e ff, log g, and perhaps abun- 
dance. Thus our initial expectation was that at low T e ff values, 
where we expect weak convective mixing (but do not usually se- 
curely detect it), all the stars might be reasonably mixed, or all 
not mixed. 

We also recall (Sec. 3.3) that the discrepancy between abun- 
dances deduced from Si n and m decreases as T e ff increases. 
Thus a possible explanation of the discrepancy we find is a 
systematic underestimate of for almost all the stars of our 



sample. To explain the discrepancies found among normal and 
HgMn stars, the values of T e ff would have to be systematically 
underestimated by an amount of the order of 1000 K, which 
would reduce the discrepancy by about 0.5 dex. For the mag- 
netic stars, the systematic error in T e ff needed would be of or- 
der 2000 K. Although such large errors are possible, they seem 
unlikely to us, and so we look for other explanations of the ob- 
served differences. 

We next consider two further possible explanations of the 
abundance discrepancy that we do detect in almost all the stars 
in our sample: non-LTE effects, and Si abundance stratification. 

5.1. Non-LTE effects in B-type stars 

Recent studies of non-LTE effects in the atmos phere of O and 
early B-type stars have been c arried out (e.g. iPrzvbilla et al.L 
12011b iNieva & PrzvbillaL I2012I) testing the predictions of LTE 
and non-LTE line synthesis codes for main sequence stars with 
T eff between about 15000 and 35000 K. They show that a hy- 
brid method, combining LTE model atmospheres with non-LTE 
line synthesis, is capable of yielding very accurate basic stel- 
lar parameters, and very satisfactory fit to most spectral lines. 
From comparisons of their non-LTE results to fully LTE syn- 
thesis, they conclude that pure LTE modelling yields meaning- 
ful results up to about a T e ff of 22000 K for a large number of 
spectral lines. However, line s with this behaviour must be care- 
fully identified. Neiner et al ] (l2012h have been able to determine 
abundances of a number of elements for the star HD 96446 (7^ 
= 21600 + 800 K) by using only spectral lines that are approx- 
imately the right strength when computed with LTE synthesis 
and non-LTE abun dances for standard stars such as a Pyx. 

Unfortunatelv. lNeiner et al. 1 (12011 find that lines of Si n and 
m mostly are not among the lines that are predicted to have the 
observed strength in standard stars by LTE synthesis. Thus it ap- 
pears that among early B type stars, non-LTE effects are proba- 
bly imp_ortantfo£Si^_These results are consistent with comments 
by IPrzvbilla et al.1 (1201 ll) who suggest that lines of silicon are 
strongly influenced by non-LTE effects above about 15000 K. 

We do not have the tools to carry out non-LTE synthesis for 
this project. However, we have included the star HD 160762 
(T e s = 17500 K) in our sample in order to explore the on- 
set of non-LTE just above the temperature range of the stars 
we are studying. For this star, INieva & Przvbillal (120121) find 
that the abundance of Si is essentially solar, with log N^/Nn = 
-4.49 + 0.05. From our LTE analysis, it is found that the values 
of Si abundance as measured with lines of Si n (-5.34) and Si in 
(-4.10) a re inconsistent with one another, and with the value 
found by INieva & Przvbillal (120121) . The fact that with suitable 
choice of lines and appropriate non-LTE synthesis computations, 
the abundance of Si, like the other elem ents in this star that were 
measured by INieva & Przvbillal (120121) . have solar abundances 
indicate that stratification is probably not significant in this star. 
Instead, the problems in this star with the lines we are using are 
probably non-LTE effects. 

We could try to test this conclusion by showing that non- 
LTE computation of our favoured lines in HD 160762 lead to an 
essentially solar, and consistent, abundance of Si. While we do 
not have the correct tools for t his computation, we ca n try to use 
the non-LTE calculations of Becker & Butler ( 1990a), who have 
computed non-LTE equivalent widths of a number of Si n and in 

1 Derived from AA4128 and 4130 simultaneously. 

2 Derived from AA5041, 5055, and 5056 simultaneously. 

3 Derived from AA4552, 4567, and 4574 simultaneously. 
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Tabl e 5. Silicon abundan ces deduced for HD 160762 from lines of Si n and m using the non-LTE equivalent width calculations of 
iBecker & Butlerl ([1990a), and LTE models computed with zeeman. 



Ion 


Line (A) 


Wa(A) 


logCAWAW 
Becker & Butler 


Zeeman 


£ = km s -1 


f = 3kms- 1 


£ = 1 km s- [ 


Sin 


4128 


57.3 


< -5.0 


< -5.0 


-5.13 1 




4130 


63.2 


< -5.0 


< -5.0 






5041 


40.1 


< -5.0 


< -5.0 


-5.34 2 




5055 


49.1 


< -5.0 


< -5.0 






5056 


21.9 


< -5.0 


< -5.0 




Sim 


4552 


70.1 


-4.1 


-4.3 


-4.10 3 




4567 


47.5 


-4.3 


-4.5 






4574 


28.2 


-4.2 


-4.5 





lines over a large grid of r e ff, log g, and We have measured the 
equivalent widths of a number of spec t ral line s in HD 160762, 
and used the tables o f lBecker & Butlerl (Q990a) to deduce corre- 
sponding Si abundances. The results of this exercise are shown 
in Table[5] together with our own LTE abundances from the same 
lines. 

Assuming the correctness of the Si abundance of -4.49, we 
find that the Si abun dance determined from li nes of Si m, using 
the computations of IBecker & Butleil (Il990al) and interpolating 
to f = 1 km s _1 , is about -4.25, somewhat closer to the correct 
value than the LTE abundance of -4.10. In contrast, the non- 
LTE abundance deduced from lines of Si n appears to be nearly 
one dex low, as is the LTE Si n value. It thus seems that the 
computations of Becker & Butler ( 1990a) for Si m multiplet (2) 
do not fully correct for non-LTE effects, and those for the optical 
lines of Si n do not even come close to correctly describing the 
non-LTE effects. (Note that the non-LTE treatment of Si n is still 
apparently unsatisfactor y today; the mode l non-L TE spectra for 
early B stars shown bv lNieva & Przvbillal (1201 2|) in their Figs 8 
and 9 do not include any of the strong optical Si n lines.) 

Our conclusion is that by an effective temperature of about 
17500 K non-LTE effects almost certainly do undermine Si 
abundance determinations from LTE spectrum synthesis, espe- 
cially for lines of Si n, but we are not at present able to quantify 
this effect, or even say how far down in effective temp erature it 
reaches. The fact that iMashonkina et al.l d2005l 120091) find im- 
portant non-LTE effects in formation of lines of Pr and Nd n and 
in in magnetic Ap stars with T e s as low as 7250 K suggests that 
non-LTE effects in Si could potentially be important well below 
the conventional limit of T e g m 15 000 K. 

5.2. Stratification 

A rather different possible explanation for the Si n - Si m abun- 
dance discrepancies is the possibility that Si is not uniformly 
abundant with height in the atmosphere of the stars studied 
here. We have already mentioned this possibility earlier when 
we pointed out that the stars without the abundance discrepancy 
all show evidence of non-zero £ values, and hence are presumed 
to have convective mixing present in the atmosphere, suggest- 
ing that these stars at least may have atmospheres that are too 
well mixed to support any kind of vertical non-uniform abun- 
dances. However, for most of the stars in this survey, we have no 
strong evidence of vertical mixing, and in the case of the mag- 
netic stars, the large energy density in the magnetic fields al- 
most certainly quenches any convective turnover, so that vertical 
abundance stratification is quite possible. 



As discussed above, vertical stratification has already 
been detected in rare e arths in cool magnetic Ap stars by 
Rvabchik ova et al.l d2001l) . who found that Pr and Nd have sub- 
stantially larger abundance relative to H in the upper atmo- 
spheres of pulsatin g magnetic (roAp ) stars. This effect has been 
further studied by iKochukhov et al.l d2009) and others cited in 
the Introduction. It now seems quite clear that stratification does 
occur to a very important extent in the cool magnetic Ap stars, 
and that it affects many chemical elements, including Si. Is there 
any evidence that stratification can explain the abundance dis- 
crepancies documented in this paper? 

We have carried out a series of simple experiments in which 
we introduce a simple smoothed step function of specified am- 
plitude and at a specified optical depth 75000 at 5000 A. This was 
done using a slightly modified version of zeeman that, given the 
optical depth and abundance change of the step, can optimise 
the fit to one or several spectral lines by adjusting the abundance 
below the step in the atmosphere. We mainly experimented on 
the two magnetic Bp stars BD+00 1659 and BD-19 5044L, both 
of which have values in the middle of our sample, around 
12 500 K, and such low values of v sin i that we can study the fit 
to line profiles, not merely to equivalent widths. This is a far sim- 
ple r scheme than the detaile d modelling described, for example, 
by IKochukhov etafl ((2006). However, it enables us to explore 
simply the plausibility of stratification as an explanation of the 
Si n-m discrepancy. 

We have explored mainly step functions in which the abun- 
dance decreases outwards. The few experiments with Si con- 
centrated high in the atmosphere lead to even larger discrep- 
ancies between the abundances determined from the two ioni- 
sation states than with an assumed uniform abundance. We have 
roughly explored abundance decreases of 1 or 2 dex, located typ- 
ically between 75000 = 0.2 and 0.8. It is found that most such 
models substantially reduce the discrepancy between the abun- 
dances deduced from the two ion states, often to below 0.5 dex. 
With a decrease of 1 dex in Si abundance at T5000 ~ 0.4 the 
discrepancy can be reduced to as low as 0.3 dex. Larger steps 
of 2 dex lead to nominal agreement between the two ionisation 
stages, but also lead to wings on the Si 11 lines that are unac- 
ceptably broad (a problem that is only visible because we are 
studying stars of low v sin i). 

The conclusion of our limited experiments is that vertically 
non-uniform Si abundance, with a substantially lower abundance 
high in the atmosphere relative to T5000 = 1> can explain at least 
a large part of the discrepancy between abundances of Si de- 
duced from the two ion states for magnetic Bp stars, and may 
entirely explain the Si n-m differences for normal and HgMn 
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stars. Clearly more detailed stratification modelling of individ- 
ual stars would be very illuminating. Similarly, further study of 
the non-LTE behaviour of the Si ions in the T e g range of 10000 
to 15 000 K would be very valuable. 

6. Discussion 

The results of our abundance analysis of the singly and doubly- 
ionised states of silicon in B-type stars are summarised as 
follows: 

1 . For a solar abundance of Si, we begin to detect Si in in main 
sequence stars at ~ 1 1000 K. 

2. All magnetic Bp stars have a discrepancy between the abun- 
dances derived from lines of Si n and Si m with the latter 
being between 0.6 -1.7 dex higher. The same discrepancy is 
observed in the non-magnetic stars, but to a much lesser ex- 
tent, with an enhancement of the Si m abundance by between 
0.3 - 0.8 dex. 

3. Only three non-magnetic stars do not exhibit a discrepancy 
between the abundance of Si derived from both ionisation 
states. These are exactly the stars where we detect a non-zero 
microturbulence parameter: HD 22136 (normal), HD 57608 
(HgMn), HD 27295 (HgMn). This suggests the presence of 
convective mixing in the atmospheres of these stars, inhibit- 
ing the development of any stratified abundance distribution. 

4. Possible errors in the atomic data (of the order of +0. 1 dex) 
cannot explain this discrepancy. Although non-LTE effects in 
Si lines may be important, it is unclear to what extent these 
will influence the derived abundances below T e g ~ 15000 K. 

5. Simple stratification models with a lower abundance high in 
the atmosphere relative to T5000 = 1 can explain a large part 
of the discrepancy between the derived abundances of Si n 
and Si m in the magnetic Bp stars, and may completely ex- 
plain the phenomena in the normal and HgMn stars. 

To fully understand the discrepancy noted between the abun- 
dances derived from the first and second ions of Si, more detailed 
stratification models are warranted. Further, studies of the non- 
LTE behaviour of Si below T e g = 15000 K would reveal to what 
extent the observed discrepancy is influenced by non-LTE ef- 
fects. 
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